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This paper studies vortex-lattice phases of rapidly rotating superfluid '^He based on the Ginzburg- 
Landau free-energy functional. To identify stable phases in the p-Q. plane {p: pressure; fi: angular 
velocity), the functional is minimized with the Landau- level expansion method using up to 3000 
Landau levels. This system can sustain various exotic vortices by either (i) shifting vortex cores 
among different components or (ii) filling in cores with components not used in the bulk. In addition, 
the phase near the upper critical angular velocity Q.c2 is neither the A nor B phases, but the polar 
state with the smallest superfluid density, as already shown by Schopohl. Thus, multiple phases 
are anticipated to exist in the p-Q. plane. Six different phases are found in the present calculation 
performed over 0.000ir2c2 < ^ < Q,c2, where Q,c2 is of order (1— T/Tc)xl0'' rad/s. It is shown that the 
double-core vortex experimentally found in the B phase originates from the polar hexagonal lattice 
near f2c2 via (i) a phase composed of interpenetrating polar and Scharnberg-Klemm sublattices; (ii) 
the A-phase mixed-twist lattice with polar cores; (iii) the normal-core lattice found in the isolated- 
vortex calculation by Ohmi, Tsuneto, and Fujita; and (iv) the A-phase-core vortex discovered in 
another isolated-vortex calculation by Salomaa and Volovik. It is predicted that the double-core 
vortex will disappear completely in the experimental p-T phase diagram to be replaced by the 
A-phase-core vortex for f2> 10'^ ~ 10* rad/s. C programs to minimize a single-component Ginzburg- 
Landau functional are available at http://phys.sci.hokudai.ac.jp/~kita/index-e.html. 



PACS numbers: 67.57.Fg, 74.60.-w 

I. INTRODUCTION 

Rotating superfluid '^He with 9 complex order parame- 
ters can sustain various exotic vortices not observable in 
superfluid ''He. This system can be a text-book case of 
unconventional vortices realized in multi-component su- 
perfluids and superconductors. I here report the richness 
and diversity of the vortex phase diagram in the unex- 
plored region of rapid rotation, wishing to stimulate ex- 
periments in the frontiers as well as to give hints to what 
may be expected in the vortex phases of multi-component 
superconductors. 

Extensive efforts have been made both theoretically 
and experimentally to clarify vortices of superfluid "^Hc; 
see Refs. 1-7 for a review. With multi-component or- 
der parameters, this system can be a rich source of un- 
conventional vortices. Those already observed in ro- 
tation up to 3 rad/s include: superfluid-core vortices 
in the B phase, i.e., the A-phase-core and double-core 
vortices;®"^^ vortices due to textures of the 1-vector in 
the A phase, i.e., the locked vortex 1,^^'^^ the continu- 
ous unlocked vortex, ^"'■^^ the singular vortex, ^^'^^ and the 
vortex sheet. The superfluid cores are possible in the B 
phase because the components not used in the bulk are 
available to flU in them. On the other hand, the A phase 
has a unique property that it can sustain vortices by a 
spatial variation of I without any amplitude reduction, as 
first pointed out by Mermin and Ho."'^ Thus, experiments 
have already revealed rich structures. 

Although SI ^ 3 rad/s is three orders of magnitude 
faster than the lower critical angular velocity fici ~ 10"'^ 
rad/s for a typical experimental cell of diameter ~5 mm, 
it is still far below the upper critical angular velocity 



^c2 ^ i^ — T/Tc) X 10'' rad/s. Thus, theoretical calcula- 
tions have been performed mostly within the isolated- 
vortex approximation in the B phase, or within 
a constant amplitude in the A phase, ^^'^^'^^'^^'^^ both 
of which are justified near fici, and not much has been 
known about the phases realized in rapid rotation. On 
the other hand, the polar state should be stable near 
ric2 at all pressures, as shown by Schopohl^^ and later 
by Scharnberg and Klemm^^ in a different context of p- 
wave superconductivity. This is because the line node of 
the polar state is most effective in reducing the kinetic 
energy dominant near f2c2- Thus, the phase near r^c2 is 
completely different from the experimentally observed A 
and B phases at $7 = 0, and there should be novel phases 
between n^i and ric2- Although Oc2 (l-T/Tc) x lO"^ 
rad/s may not be attainable in the near future, clarifying 
the whole phase diagram of f2ci ^ ^ ^c2 would stim- 
ulate experimental efforts towards the direction; it cer- 
tainly remains as an intellectual challenge. In addition, 
such a study will be useful in providing an insight into 
the vortices of multi-component superconductors where 
Hc2 can be reached easily. 

Following a previous work,^^ I present a more exten- 
sive study on vortices in superfluid '^He. To this end, I 
adopt the standard Ginzburg-Landau (GL) free-energy 
functional valid near Tc, as most calculations performed 
so far. To clarify the vortex phase diagram, however, I 
take an alternative approach to start from r2c2 proceed- 
ing down towards fJci as close as possible. A powerful 
way to carry out this program is the Landau-level ex- 
pansion method (LLX), developed recently^^'^^ and ap- 
plied successfully to several other systems. Combin- 
ing the obtained results with those around f2ci will pro- 
vide a rough idea about the whole phase diagram over 
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rJci <^< ^c2- It should also be noted that the results 
from the GL functional are expected to provide qualita- 
tively correct results over <T <Tc, as supported by a 
recent isolated- vortex calculation on the B phase^° using 
the quasiclassical theory.^^ 

This paper is organized as follows. Section II presents 
the GL functional. Section III explains LLX to mini- 
mize the GL functional. Section IV provides a group- 
theoretical consideration on the classification of vortex 
lattices and the phase transitions between them. Sec- 
tions V presents the obtained p-0, phase diagram over 
0.01f7c2 <^< f^c2 together with detailed explanations on 
the phases appearing in it. Section VI discusses a phase 
transition between the A-phase-core and double-core lat- 
tices extending the calculation down to 0.0001Qc2- Sec- 
tion VII summarizes the paper. Appendix A presents 
basic properties of the basis functions used in LLX. 



II. GINZBURG-LANDAU FUNCTIONAL 

Superfluid ^He is characterized by 9 complex order pa- 
rameters A^i (/i, i = x,y, z) inherent in the p-wavc pairing 
(X = 1) with spin 5'= 1, where /x and i denotes rectangular 
coordinates of the spin and orbital spaces, respectively. 
The GL free-energy functional near Tc is given with re- 
spect to the second- and fourth-order terms of A^j. Using 
the notation of Fetter/ the bulk energy density reads 

'^f^'i^t^i^vi^M^id +/^4A^^Aj^^-A^jA^i 
'^Ps^fj.i-^fj.j^ui^i^j ) (1) 

where a and ,6j arc coefficients, and summations over 
repeated indices are implied. The gradient energy density 
is well approximated using a single coefficient K as 

/k = K{d:A;,djA^j+d:A;jdiA^j+d*A;jdjA^i) , (2) 

where d is defined in terms of the angular velocity as 

.2m3 



d = V 



'-flxr. 



(3) 



In addition, there are tiny contributions from the dipole 
and Zeeman energies: 

/m = QmllfiA^^HiiAiii , (5) 

respectively. Given the coefficients in Eqs. (l)-(5), the 
stable state can be found by minimizing 



F = Fo+Fi 



(6) 



with 



Fo = ^J{h+h)dv, (7a) 
F^ = ^liU+fm)dr, (7b) 



where V is the volume of the system. Important quanti- 
ties in the functional are 



e^(x/a)i/2 = e(0)/(i-r/r,)i/2 



(8b) 
(8c) 



which define the GL coherence length, the dipole length, 
and the characteristic magnetic dipole field, respectively. 

The coefficients a, (3j, K, g„^, and are fixed by ex- 
actly following Thuneberg's procedure^^ used in identi- 
fying the B-phase core structures as follows. The weak- 
coupling theory yields^ 



. = «(i-r/rj, 



w . 



(9a) 



120(7r/cBTc)2 



K = 



240(7rfcB7;)2 



(9c) 



where A^(0) and vp are the density of states per spin 
and the Fermi velocity, respectively. The coefficients a 
and K are estimated by Eqs. (9b) and (9c) using the 
values of A''(0), Tc, and vp determined experimentally 
by Grcywall.'^^ In contrast, cannot account for the 
stability of the A phase at f2 = 0. Strong-coupling cor- 
rections are included in (3j by (i) using the Sauls-Serene 
values for 1.2 MPa<p<3.44 MPar"^^^ (ii) adopting the 
weak-coupling result at p = MPa; and (iii) inter- 
polating the region MPa < p < 1.2 MPa. With this 
procedure, the A-B transition for il = is predicted at 



Ppcp ■ 



= 2.85MPa. 



(10) 



It thus yields a qualitatively correct result that the A 
phase is stabilized on the high-pressure side, although 
the value is slightly higher than the measured 2.1 MPa. 
The values of ga have been studied extensively in a recent 
paper by Thuneberg.^^ It is shown that the following for- 
mula nicely reproduces the values extracted from various 
experiments: 



Mo 

9'^ = ^ 



jhN{0) In 



1.1339 x0.45Tf 



(11a) 



where /jq and 7 denotes the permeability of vacuum and 
the gyromagnetic ratio, respectively, and Tp is the Fermi 
temperature defined by Tp = 3n/4A^(0)fcB with n denot- 
ing the density. As for g^, the following weak-coupling 
expression is sufficient: 



7a:y)N ((})(' 
48[{1+ FS)7rkBTc] 



(lib) 



with Fq the Landau parameter. The values of F^ are 
taken from Wheatley'^^ but corrected for the newly de- 
termined Ff by Greywall;^^ this F^{p) is tabulated con- 
veniently in Rcf. 4. Table I summarizes the pressure de- 
pendences of basic quantities thus calculated. 
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TABLE I: Numerical values for the parameters of the GL theory at different pressures. See text for details of the calculations. 



p 


a/il-T/Tc) 
[10^0 J-im-3] 




K 


£(l-T/TcY''^ 




Ha 


fic2/(l-T/rc) 


[10® Pa] 


[lO^s J-3m-3] 


[10^4 J-^m-i] 


[10-® m] 


[10"^ m] 


L ^ J 


[10® rad/s] 


0.0 


1.68 


21.7 


41.8 


5.00 


4.91 


0.803 


2.10 


0.4 


2.03 


11.5 


17.8 


2.96 


2.85 


1.09 


5.98 


0.8 


2.33 


8.55 


11.2 


2.19 


2.05 


1.32 


10.9 


1.2 


2.59 


7.29 


8.32 


1.79 


1.64 


1.48 


16.4 


1.6 


2.85 


6.71 


6.71 


1.54 


1.38 


1.60 


22.3 


2.0 


3.09 


6.43 


5.70 


1.36 


1.19 


1.73 


28.5 


2.4 


3.33 


6.33 


5.02 


1.23 


1.05 


1.85 


34.9 


2.8 


3.57 


6.39 


4.54 


1.13 


0.949 


1.96 


41.3 


3.2 


3.81 


6.56 


4.19 


1.05 


0.864 


2.04 


47.8 


3.44 


3.96 


6.69 


4.02 


1.01 


0.821 


2.10 


51.7 



It should be noted finally that, although Thuncberg's 
procedure is adopted here, qualitative features of the 
phase diagram (Fig. 1) obtained below are expected to be 
the same among the models for f3j = j3j (p) which yields 
the A-B transition at f2 = 0. This has been checked 
for the spin-fluctuation-feedback model of Anderson and 
Brinkman.^^'^^ 



III. METHOD 

To minimize Eq. (6), let us first rewrite the order pa- 
rameters as 



-R A ■ 



(12) 



where R^i, denotes the spin-space rotation and A^i is the 
order parameter of the restricted space where the spin 
coordinates arc fixed conveniently relative to the orbital 
ones. We then find that the gradients of A^i and R^v 
do not couple in Eq. (2) due to the orthogonormality: 
RfixRi^x = ^ij,u- The characteristic lengths for Ai,i and 
R^i, are given by ^ and of Table I, respectively, where 
we observe that is much longer than both ^ and the 
intervortex distance 1^ defined below in Eq. (14) for the 
relevant range 0.0001f2c2 < ^ < ^c2 considered here. It 
hence follows that R^^ is virtually kept constant in space. 
It can be fixed by Eq. (7b) after obtaining from Eq. 
(7a). This is due to the smallness of Eq. (7b) relative to 
Eq. (7a). 

To minimize Eq. (7a) with respect to I assume 

uniformity along $7 || z. I then define creation and anni- 
hilation operators {a) ,a) satisfying aa^ — a)a=l as 



^/2 



V2 



(d^+tdy), (13) 



with 



(14) 



It is also convenient to introduce the quantities: 



iW^V, 4±i)^i=(4.Tii^,), (15) 



which denote the expansion coefficients of and {kxi 
iky) 1^/2, respectively. Now, Eq. (2) can be rewritten in 
terms of Eqs. (13)-(15) as 

/k= |{(l + |m|) [i(")*i(r^+2(ai(r^)*ai(r)] 

-2[(aiW)*ati(-i)+„i(i)(„ti(-i))*]}. (16) 
Equations (1), (4), and (5) arc transformed similarly v&- 



ing A^jAj/j = A\, ' and A^.iAi,i= A'li, 'A\, 

Prom Eq. (16) and the bilinear term in Eq. (1), we find 
that the p-wave superfluid transition is split, due to the 
uniaxial anisotropy Ct, essentially into those of different 
m channels as 



(m) Ai-m) 



"c2 ^ 



f2 



(-1) 

c2 



ha 
AmzK 

3+V6 



c2 



c2 



(1) 



1 



c2 



polar 

ABM-(SK) 
ABM+ 



(17) 



in agreement with the results of Schopohl^^ and Scharn- 
berg and Klemm.^* Here the first and third ones cor- 
respond to the polar state (oc k^) and the Andcrson- 
Brinkman-Morel (ABM) state (oc kx + iky) with 1^ = 1, 
respectively, both in the N = Landau level. In con- 
trast, the ABM~ (also called as Scharnberg-Klemm or 
SK state) is not a pure ABM state; it is made up of the 
kx — iky state in the N = Landau level and the kx+iky 
state in the N = 2 Landau level. Equation (17) tells us 

that it is the polar state A^^^ j^O which is realized at i}c2- 
Its stability is due to the line node in the xy plane perpen- 
dicular to ft which works favorably to reduce the kinetic 
energy dominant near Qc2- Indeed, given the condition 
that the average gap amplitude on the Fermi surface is 
constant, the superfiuid density pxx of the polar state 
oc3 / dOkfc^A:^ = 3/15 is half that of the ABM state 
2 / dQ]ikl{kl + ky) — 6/15. This explains the 

difference of the factor 2 between Oc2 — f^c? ^^'^ ^c2^- 
As seen in Table I, this Qc2 is of order (1-T/Tc) x lO'^ 
rad/s. 



„polar 
fxx 
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With these preparations, Eq. (7a) is minimized by 
LLX.^® In the end, it can be performed quite efficiently, 
especially near flc2, with the Ritz variational method of 
expanding A^™' in some basis functions and carrying out 
the minimization with respect to the expansion coeffi- 
cients. Convenient basis functions for periodic vortex 
lattices are obtained using the magnetic translation op- 
erator: 



as 



rR = cxp[-R- (V + 2im3nxr/h)] , 



(18) 



where R denotes a lattice point spaned by two basic vec- 
tors: 



ai = (aix,aiy,0) 
a2 = (0,02,0) 



01x02 = 27rL . 



(19) 



Hence those lattices have a unit circulation quantum 

k = /i/2to3 per every basic ccU,^ as required. The basic 
vectors of the corresponding reciprocal lattice are given 

by 



bi = (a2Xz)//2 
b2 = (zxai)/Z2 



(20) 



The magnetic Bloch vector is defined in the Brillouin 
zone of the reciprocal lattice as 



q= -r7bi + —rh2 



(21) 



where Aff is an even integer with Mf denoting the number 
of K in the system. Using these quantities, the basis 
functions are obtained as simultaneous eigenstates of a^a 
and Tr as^^ 



E 

n=-A^f/2-|-l 



xe 



2TTlc/a2 (x-llqy-naix 



H 



N 



Ic. 



(22) 



where N denotes the Landau level and Hn is the Hcrmitc 
polynomial.^* Useful properties of V'JVq are summarized 
in Appendix A. Let us expand ^IT"^ in {f/'Afq} as 



i(r)(r)=x/y 5:5:47)«r). 



(23) 



N=0 q 



Then substitute Eq. (23) into Eq. (7a), perform a change 
of variables {x,y) = {aixS,a\yS + a2t), and carry out the 
integration in Eq. (7a) with respect to (s.t). Now, Eq. 
(7a) is transformed into a functional of the expansion 
coefficients {c|^™^} and a couple of lattice parameters: 



f3 = cos ^{aiy/ai) : apex angle 
p = 02/01 : length ratio 



(24) 



(25) 



For a given fl, this Fq is minimized directly. 

Numerical minimizations have been performed as fol- 
lows: (i) Cut the series in Eq. (23) at some A^max and sub- 
stitute it into Eq. (7a). The convergence can be checked 
by increasing A^,„ax, as increasing A^max is guaranteed to 
yield a bc;tter solution with a lower Fq. (ii) Numerical 
integrations in Eq. (7a) are performed using the trape- 
zoidal rule which is known very powerful for periodic 
functions. To this end, the basis functions of A^< A^max 
on the discrete points are tabulated at the beginning of 
the calculations. Equation (25) is thereby evaluated for 
each set of arguments. It has been necessary to increase 
Aniax from 6 near 0^2 to 360 at 0.01flc2, and the inte- 
gration points in the basic cell from 6^ to 36^ accord- 
ingly to obtain the relative accuracy of order 10~* for 
Fq. (iii) Equation (25) is minimized first with respect 

to {c^^'} by Powell's method or the conjugate gradient 

method, and then, if necessary, with respect to (/?, p) 
by Powell's method. The conjugate gradient method is 
about 10 times faster, although the programming is far 
more cumbersome. Powell's method is fast enough for 
ri>0.05ric2- (iv) Second-order transitions are identified 
carefully by combining group-theoretical considerations 
of Sec. IV with signals of the relative change in the slope 
of dFo /dfl. Its analytic expression is obtained as^^'*° 



^ I f ^ 



(26) 



so that it can be calculated quite accurately without re- 
course to any numerical differentiation. 

Once Af^i are fixed as above, Eq. (12) is substituted 
into Eq. (7b). Considering i?^^ constant in space, as 
mentioned before, the three independent parameters of 
the matrix R are then determined by minimizing Eq. 
(7b). 

Some phases encountered below have a common fea- 
ture that all the components other than 

^={,RxZ1 RyZI Rzz) \ A= (Aj;x) ^22) (27) 

can be put equal to zero. For those states, Fx is given by 
Fx = gd{d^M^J,-l{A*-A)) , (28) 
where {■■ ■) = y J ■ ■■ dr, and M^y is defined by 

M^, = {2A;A, + k-AH^H,/Hl) . (29) 

It hence follows that d points parallel to the eigenvector 
corresponding to the smallest eigenvalue of M. 
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IV. GROUP-THEORETICAL 
CONSIDERATIONS 

A. Classification of Vortices 

As in the case of ordinary solids,^^'^^ various vortex 
lattices can be classified according to their symmetry. It 
turns out that the operator (18), the basis functions (22), 
and the expansion (23) are quite useful for this purpose. 

Classification of isolated vortices in superfluid ^He 
has been carried out by Salomaa and Volovik^'^° and 
Thuneberg.-^^ Such classification for vortex lattices has 
been performed recently by Karimaki and Thuneberg,^^ 
without recourse to Eqs. (18)-(23), however. It is shown 
below that making use of Eqs. (18)-(23) provides a trans- 
parent classification scheme. 

To start with it is necessary to define "symmetry op- 
erations" unambiguously, since we arc considering a phe- 
nomenological GL functional rather than a microscopic 
Hamiltonian. In the case of a Hamiltonian, "symmetry 
operations" mean those operations which commute with 
the Hamiltonian. In the present case of a GL functional, 
"symmetry operations" are defined as those operations 
which keep all the physical (i.e. observable) quantities 
invariant. Let us restrict ourselves to orbital degrees of 
freedom, as appropriate when Eq. (7b) can be regarded 
as a tiny perturbation. Then, besides /b+/k of Eqs. (1) 
and (2), there are three relevant physical quantities: 



p(r) E 



A* A ■ 



(30a) 



h 



im(A*,aA 



+Al^d,A^i) , (30b) 

h{r) = -ieijkA*^jAf,k/p{r) , (30c) 

which denote the density, current, and normalized or- 
bital angular momentum of Cooper pairs, respectively. A 
general operator <S transforms A^^j, di, and Sijk in these 
expressions as 



di - 



9. 



(5) 



^iifc = Sii'Sjj'Skk'^i'j'k' , 



(31a) 
(31b) 
(31c) 



where (5)^,^ = 5^,^ denotes a matrix representation of 
S. Now, symmetry of a vortex lattice is specified by 
those operations {<S} which keep the physical quantities 
invariant as 



/b(r) + /k(r) ^ 

j(r)=j('5Hr) 
I(r) =l('S)(r) 



/r'w+/f^(r) 



(32) 



where j'-'^-', for example, denotes the expression obtained 
by substituting Eq. (31) into Eq. (30b). The collection of 



these operations {S} forms a group, which characterize 
the relevant vortex lattice. 

Since Tr of Eq. (18) commutes with d of Eq. (3), the 
operator Tr plays exactly the same role as the transla- 
tion operator of the ordinary space group. This is one 
of the main reasons why the basis functions (22), which 
are eigenstates of Tr, is advantageous in expanding the 
order parameters as Eq. (23). 

As an example, let us specifically consider Phase IV 
of Fig. 1 given below. It corresponds to the hexagonal 
lattice of the normal-core vortex in the B phase which was 
discovered in the isolated-vortex calculation by Ohmi, 
Tsuneto, and Fujita.^^ The present calculation with LLX 
shows that non-zero order parameters can be expressed 
as 



(33a) 



N 



4"'^ = ^ Et4q>iVqi +4+iq,^iV+4q.] , (33b) 



JV 



= iW 5^[4q>iVqi -4+U,^iV+4q,] , (33c 



N 



4'^ = Et WiVqi+4V2q,^^+2q,] , (33d) 



JV 

= -iW E[4U^qi -4V2q,V'iV+2qJ , (33e) 
N 

where the summations run over N = 6n (n = 0, 1, 2, • • • ), 
qi is a magnetic Bloch vector which can be chosen ar- 
bitrarily due to the broken translational symmetry, and 
all c^'^'s are real. It is convenient to fix qi = i(bi-|-b2) 

so that a core is located at the origin, i.e., if"^(O) = 0. 
Then, with the properties of i^Nqi given by Eq. (A8), one 
can show that the order-parameter matrix (A)^j = A^i 
satisfies 



R^A{B^'r)R^'=e-"'A{r), 
^.iV^^r)^-i=e-5-/*i(r), 



^j;iV/r)a/ =e- 



"/^i(r 



(34a) 
(34b) 
(34c) 
(34d) 



where and {i = x, y) denote, respectively, a rotation 
around the z axis by <y9 = n7r/3 and the mirror reflection 
with respect to the plane perpendicular to the i axis. Us- 
ing Eq. (34) in Eq. (31), one can show that Eq. (32) holds 
for the four operations of Eq. (34). Similar considerations 
show that, besides Tr, there are symmetry operations in 
Phase IV listed in the fourth row of Table II. Here the 
primitive cell is spanned by the basic vectors (19), thus 
containing a single circulation quantum. If Tr is identi- 
fied with the usual translation operator, this space group 
can be labeled as PG/mw'm'.^^ 

The same consideration is performed for every stable 
vortex phase found numerically. The results are summa- 
rized below in Table II. 
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B. Phase Transitions 

The expansion (23) is also useful for enumerating pos- 
sible transitions in vortcx-latticc pliascs. Indeed, we 
can use the known results from the space group^^ us- 
ing the correspondence exp(— R- V) ^ Tr and (Bloch 
states) ^(magnetic Bloch states). 

We first summarize basic features of the conventional 
Abrikosov lattices with a single order parameter within 
the framework of LLX:^^ (a) Any single q = qi suffices 
to describe them, due to broken translational symmetry 
of the vortex lattice. Each unit cell has a single circu- 
lation quantum k. (b) The hexagonal (square) lattice is 
made up of A'' = 6n (4n) Landau levels (n = 0, 1, 2 • • • ), 
and the expansion coefficients c^cn can be chosen real, 
(c) More general structures are described by A^ = 2n lev- 
els. Odd-A' basis functions, some of which have finite 
amplitudes at the cores of even-A basis functions (see 
Appendix A), never participate in forming the order pa- 
rameter, because such mixing is energetically unfavor- 
able. 

With multi-component order parameters, there can be 
a wide variety of vortices, which may be divided into two 
categories. I call the first category as "fill-core" states 
with a single circulation quantum per unit cell, i.e., only 
a single q = qi is relevant. Here the cores of the conven- 
tional Abrikosov lattice are filled in by some superfluid 
components using odd- A' wavefunctions of Eq. (22). The 
second category may be called "shift-core" states, where 
core locations are not identical among different order- 
parameter components. The corresponding lattice has 
an enlarged unit cell with multiple circulation quanta. 
General structures with circulation quanta per unit 
cell can be described using different q"s. and the unit 
cell becomes times as large as that of the conven- 
tional Abrikosov lattice. For example, structures with 
two quanta per unit cell can be described by choosing 



qi 



bi + b2 



q2=o, 



(35) 



where bi and b2 are reciprocal lattice vectors defined by 
Eq. (20). 

With these observations, we now realize that following 
second-order transitions are possible in a two-component 
system, (i) Deformation of a hexagonal or square lat- 
tice within a single-q lattice. This accompanies entry of 
new even-A^ basis functions, and the coefficients {cNq} 
become intrinsically complex, (ii) The entry of odd-A^ 
basis functions within a single-q lattice, i.e. a transition 
into a fill-core lattice. Here, cores of even-A' basis func- 
tions are filled in by some superfluid components not used 
in the bulk phase. The transition occurs below some crit- 
ical angular velocity smaller than f2c2/3. (iii) Mixing of 
another wave number q2, i.e. a transition into a shift- 
core lattice. When applied to this system, the Lifshitz 
condition^^ concludes exclusively that q2 — qi should be 
half a basic vector b^ (j = 1, 2) of the reciprocal lattice, 
i.e. the only possibility is doubling of the unit cell. 



Though not complete for the present nine-component 
system, the above list covers most of the transitions found 
below. The case (ii) corresponds to superfluid-core states 
such as the A-phase-core and double-core vortices in the 
B phase, whereas (iii) will be shown to describe a unit cell 
with two circulation quanta like the continuous unlocked 
vortex of the A phase. ^ Since no odd-order couplings exist 

between and (ii"^\ il,^') in Eqs. (l)-(5), we also 

realize that a state of Ai~^^ ^0 or A^^^ ^0 accompanies 
a second-order transition as we decrease ft from f2c2- 



V. PHASE DIAGRAM FOR 0.0inc2<n<fic2 

Figure 1 displays the obtained p-Q phase diagram for 
O.Oir2c2 < < ^c2, where five different phases I-V are 
present. Symmetry properties of each phase, which have 
been clarified by group-theoretical considerations similar 
to that given around Eq. (34), are summarized in Table 
n. Also, critical values of the transitions are listed in 
Table HI. Each phase is explained below in detail. 

According to Eq. (17), the normals superfluid transi- 
tion in rotation occurs at i}c2 — ^lS into Jii^^ ^0. Thus, 
in Phase I just below fic2, the polar state is realized to 
form the conventional hexagonal lattice as 



(36) 



N 



with N = 6n (n = 0, 1, 2, • • • ), /3 = 7r/3, and p = I. All 
"-TVcji ^® P^* veal, and qi is conveniently chosen as 
Eq. (35) so that a core is located at the origin. As already 
discussed below Eq. (17), the stability of the polar state 
near flc2 can be attributed to its line node in the plane 
perpendicular to which works favorably to reduce the 
kinetic energy dominant near flc2- It follows from Eq. 
(28) that the d- vector of Eq. (27) at H = lies along an 
arbitrary direction in the xy plane which is fixed sponta- 
neously. This degeneracy is lifted by a field in the plane 
so that d_LH is realized. 
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FIG. 1: The p-fi phase diagram for 0.01 < n/fic2 < 1.0. There 
are five different phases I-V. The III- IV transition is first-order 
corresponding to the A-B boundary in finite f2, whereas the 
others are second-order. 
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TABLE II: Properties of Phases I-V. The notation of Ref. 42 is used for the point-group operations with i = 1,2, 3, m = x,y, 
and s = a,b, and 6 denotes the time-reversal operation. The point-group operations of Phases I and IV- VI are defined with 
respect to the origin. In contrast, those of Phases II and III are about the point (ai+a2)/2 where a core of the SK sublattice 

is located, with x and y axes taken along ai+a2 and ai— a2, respectively. The operator {/|ai}, for example, denotes successive 
operations of the inversion and a magnetic translation by ai, i.e., {J|ai} = Tai/. 

p symmetry operations 

1 E, Cg*\ C^fK C2z, I, Sgf-*, S'g*-', ah, OC'^i, OC^'i, OcTdi, Offvi- 

1 E,C2z, {/|ai}, {cThlai}, e{C2m\&i},eam. 

1 E,C[^\C2., {/|ai}, {5<±'|ai}, {a;,|ai}, e{C2^|ai}, e{C23|ai}, ^a^, Oaas- 

1 E, Cg^', Cg^', C2, /, S*^^ , ah, 0C'2i, 0C2i, 6(Tdi, Oavi- 

1 E, Cg^\ C'aJ'; C2z, Offdi, OcTvi- 
~1 E, C2z, Oadm- 



Phase 


unit cell 


/? 


I 


ai, a2 


7r/3 


II 


ai+a2, ai— a2 


l<(3< 


III 


ai+a2, ai— a2 


it/2 


IV 


ai, a2 


tt/S 


V 


ai, a2 


7r/3 


VI 


ai, a2 


-7r/3 



TABLE III: Critical angular velocities of the phase transitions in units of Q.c2- 

p (MPa) 04 08 L2 1^6 2^0 2.4 2.85 3.2 3.44 

m 0.8527 08548 08569 08588 O8605 (18624 08644 0.8671 0.8696 0.8716 

II- III 0.8388 0.8411 0.8433 0.8454 0.8473 0.8493 0.8515 0.8544 0.8572 0.8594 

III- IV 1.18 X IQ-^ 9.00 X 10"^ 6.54 x 10"^ 4.56 x 10"^ 3.16 x 10"^ 1.89 x 10"^ 8.32 x 10"^ 0.0 — — 

IV- V 2.23 X 10"^ 1.78 x 10"^ 1.44 x 10"^ 1.20 x 10"^ 1.05 x 10"^ 9.08 x 10"^ 7.78 x 10"^ 0.0 — — 



In Phase II, ^i^^^ also become finite as^^ 

^i"'^ = ^^ E 4"qi i^N^. , (37a) 

N 

N 

The TV's in Eqs. (36)-(37) are even numbers, and all the 

coefficients are essentially complex except Cgq^ which can 

be chosen real using a gauge transformation. The I-II 
transition is second-order corresponding to the SK line 
of Eq. (17). Indeed, A^ = and TV = 2 Landau levels are 
dominant in Eqs. (37a) and (37b), respectively. Due the 

presence of A^z \ however, the critical angular velocity 
q(i<->ii) somewhat lowered from = 0.9082f2c2 of 

Eq. (17), as seen in Fig. 1. The vector q2 is shifted from 
qi by half a basic vector bj , so that the unit cell is dou- 
bled carrying two circulation quanta. Thus, Phase II is 




0.850 0.855 0.860 0.865 



FIG. 2: The apex angle /? as a function of f2/ftc2 at p = 2.0 
MPa. Phase transitions occur at Q^''"''' = 0.8623f2c2 and 
n("*"™)=0.8493nc2. 



composed of interpenetrating polar and SK sublattices. 
This superposition with shifted core locations is energet- 
ically more favorable than that with the identical core lo- 
cations, because A*^^A^i becomes finite everywhere. It is 
convenient to choose q2 = 0.''^ Then p=l holds through- 
out, and as seen in Fig. 2 calculated for p=2.0 MPa, the 
apex angle (3 changes continuously throughout the phase 
between 7r/3 and 7r/2. This is a centered rectangular lat- 
tice with the primitive vectors given by ai±a2. Figure 3 
displays the order-parameter amplitudes at p = 2.Q MPa 
over —02 <x,y <a2 for (a) fl = 0.863Oc2 with (3 = n/S 



0.0 0.06 0.12 




(a) £2=0.863Qc2 (b) i>=0.852£2g2 (c) £2=0.849Qc2 

^l!li(r>4Lu(r) at p=2.0MPa 



FIG. 3: The amplitude A'^^Af^i normalized by the f2 = value 
over -02 <x,y <a2 at p = 2.0 MPa. (a) At = 0.863nc2 
where (3 = Tr/3 (polar state), (b) At f2 = 0.852f2c2 where 
= OAUtt (Phase II). (c) At = 0.8490^2 where (3 = -k/2 
(Phase III). Below O = 0.86230^2, the SK sublattice grows 
gradually with the largest amplitude at the cores of the polar 
sublattice. Thus, the amplitude A'l^iA^i is finite everywhere 
in (b) and (c). The apex angle changes continuously from 
/9=7r/3 at = 0.8623nc2 to /3 = 7r/2 at O = 0.8493nc2. 
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phase III (p = 2.0 MPa; £1 = 0.05 £1^^) 
0.0 0.5 1.0 -1.0 -0.5 





(a)A;,(r)A^,(r) 



(b)/,(r) 




(c)/ir) 



(d)(/ir),/(r)) 



FIG. 4: Characteristic quantities in Phase III calculated over 
-02 < X, y < a2 for = 0.05r2c2 at p = 2.0 MPa. (a) A'^.A^i 
normalized by the 51 = value, (b) /z(r). (c) /a,(r). (d) 
Projection of I(r) onto the xy plane. The plot for ly{r) is 
obtained by rotating Fig. 4(c) by — 7r/2. The p and x sites 
in Fig. 4(d) correspond to the pure polar state and the pure 
ABM state with 1^ = —!, respectively. 



(still in polar state), (b) ^ = 0.852^1^2 with ^ = 0.4147r 
(Phase II), and (c) $7 = 0.849^02 with /3 = 7r/2 (just be- 
low the II-III boundary). We observe that the SK state 
grows gradually from Q = 0.8623rJc2 with the largest am- 
plitude at the cores of the polar state. The d-vector has 
the same character as Phase I. 

In Phase III, the system remains in the square lattice 
with (3 = n/2 and p=l. Only N = 4n, An, and 4n-|-2 
(n = 0, 1, 2, • • • ) Landau levels are relevant in Eqs. (36), 
(37a), and (37b), respectively. Whereas all c^q^'s can be 
put real, the coefhcients c^^^l are complex with a com- 
mon phase — 7r/4 relative to Cpq^. The II-III transition 
is second-order, as may be realized from the square-root 
behavior of Fig. 2. As is decreased from the bound- 
ary, the SK sublattice with the coefficients c^^^l grows 
rapidly. Also within the SK sublattice, the ABM(+) com- 
ponent with the coefficients c^^^ becomes less important 
at lower 17 so that the sublattice approaches to the pure 
ABM state with 1^ = —1. Figures 4(a)-(d) display the 
order-parameter amplitude (30a) and the orbital angu- 
lar momentum (30c) calculated over —02 <x,y<a2 for 




phase IV 
p = 0.0 MPa 
Q = 0.025Qc2 



ImA 



'xy 



5 

X/S, 



10 



FIG. 5; Order parameters of Phase IV along x axis at p = 0.0 
MPa and = 0.025f2c2 • The amplitudes are normalized by the 
f2 = value of A^x, and the common phase tt/8 is subtracted 
from Ai^i. 



Q = 0.05fic2 at p = 2.0 MPa. We realize from these fig- 
ures that Phase III at lower Q is essentially the A-phase 
mixed-twist lattice with polar cores and double circula- 
tion quanta per unit cell.^'^'' Indeed, as seen in Fig. 4(d), 
the 1-vector rotates from downward at the origin to hori- 
zontally outward or inward towards polar cores. A group- 
theoretical consideration, similar to that given around 
Eq. (34), clarifies that this phase is characterized by the 
symmetry operations given in the third row of Table II; 
they are defined with respect to a polar core, i.e., p site 
of Fig. 4(d). It is worth comparing Fig. 4 with Fig. 3(c) 
at f2 = 0.849ric2- The cores at 17 = 0.849^02 now acquire 
a large ABM amplitude with = — 1, and the points with 
the maximum polar amplitude at 17 = 0.84917c2 turn into 
cores, i.e., singular points of the vector field l(r) where 
the amplitude A'^^Af^i is also smallest. Notice that this 
phase is stable at all pressures. It is also worth pointing 
out that, although both carrying double quanta per unit 
cell, this structure is completely different from the vortex- 
sheet-like structure found in a two-component system. 
The d-vector has the same character as Phase I. 

Phase IV is a hexagonal lattice of the B-phase normal- 
core vortex or o-vortex, which was found in the isolated- 
vortex calculation by Ohmi, Tsuneto, and Fujita.^^ How- 
ever, it was concluded later by another isolated- vortex 
calculation^^ that this normal-core vortex is a metastable 
state. Indeed, the vortex has never been observed exper- 
imentally. The present calculation shows, however, that 
it will be stabilized as we increase 17. The non-zero com- 
ponents besides Eq. (36) are given by 

4"'^ = l][4"q>iVq, +4-+lq,V'iV+4qJ , (38a) 

JV 

= iW Yli&N^^ -4"-;lq,^iV+4q,] , (38b) 



JV 



-'Af-|-2qi 



qiJ ' 



(38c) 



JV 



iW = -iW ^[cW ,^^JVq, -4;2q,V'JV+2qJ , (38d) 



JV 
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FIG. 6: Rotation angle 9 defined by Eq. (39) as a function of 
ft/f2c2- It is calculated beyond the III- IV boundary in Fig. 1 
into the metastable region. 

where A'' = 6n, /3 = it/3, p= 1, and all the coefficients 
arc real. Figure 5 plots the order-parameter amplitudes 
along X direction calculated for il ~ 0.025r2c2 at p = 0.0 
MPa. They already display the same features as those 
obtained by the isolated-vortex calculations. The 
III-IV boundary in Fig. 1 is a first-order-transition line, 
approaching towards p = Ppcp at = and vanishing for 
P>Ppcp, as expected. Thus, this line may be regarded 
as the A-B phase boundary in finite O. It is convenient 
here to parameterize the rotation matrix R of Eq. (12) 
as^ 

Rnv = S^v COS 6 +n^nu{l- COS 9) - s ^vxn\ sm9 , (39) 

where n is a unit vector. Then it is found that, as 
is decreased, the rotation angle 9 approaches from below 
to the bulk value cos~^(— 1/4), as seen in Fig. 6. The n- 
vector at iJ = lies along an arbitrary direction in the xy 
plane which is fixed spontaneously. Upon applying the 
magnetic field in the plane, the state n|| H is realized. 

Phase V continuously fills in Phasc-IV core regions 
with the superfiuid components not used in Phase IV. 
Besides those of Phase IV, the solution indeed has new 
non-zero components which are made up only of odd Lan- 
dau levels as 

= ^ El ^i^VlqiV'iV+lqi +ci^l5qiV'iV+5q.] , (40a) 
N 

if = iW E[c^llqiV'iV+lq.-C^V5q,V'iV+5qJ , (40b) 

N 

4"'^ = ^ E 4"+ki^^+5qi ' (40c) 

JV 

= E ci^+lqxV'iV+lq, , (40d) 

JV 

where N = 6n, /3 = 7r/3, p = 1, and all the coeffi- 
cients are real. Figure 7 plots the order-parameter am- 
plitudes along X direction calculated for O = 0.0ir2c2 at 
p = 0.0 MPa. The core superfiuid components satisfy 
KeA^x = ImAzy and KeA^z = ImAy^ at the origin. As 



FIG. 7: Order parameters of Phase V along x axis at p = 0.0 
MPa and n = 0.01f2c2- The amplitudes are normalized by the 
n = value of Axx, and the common phase 7r/8 is subtracted 
from A^i. The equalities ReA^^! = Im^^j, and Reji^^ = Im^yj 
hold at x = 0. 



may be realized from this figure, this state corresponds to 
the A-phase-corc vortex or axisymmctric v-vortex foimd 
theoretically by Salomaa and Volovik,^° which also has 
been observed experimentally.^'^ Indeed, using the prop- 
erties of ipNqi given in Eq. (A8), one can show that the 
order-parameter matrix A = (A^i) satisfies Eqs. (34a)- 
(34c) but not Eq. (34d). Thus, the present calcula- 
tion has clarified for the first time that a vortex lattice 
with superfiuid cores can be described by a superposi- 
tion of odd Landau levels. The IV-V phase boundary 
is a second-order transition, as anticipated by Salomaa 
and Volovik based on a single- vortex consideration,^*^ 
which is driven mainly by the N = 1 Landau level. Due 
to the presence of finite order parameters composed of 
even Landau levels, however, the critical angular veloc- 
ity Q^c^'^^'' is lowered from 0^2/3 expected for the pure 
N = l Landau level to become smaller than ~ 2xl0~^f2c2 
rad/s, as seen in Fig. 1. The n- vector has the same char- 
acter as Phase IV. 



VI. TRANSITION BETWEEN A-PHASE-CORE 
AND DOUBLE-CORE LATTICES 

Using the same parameters as those given in Table I, 
Thuneberg^^ carried out an isolated-vortex calculation. 
He thereby succeeded in identifying two kinds of vortices 
experimentally found in the B phase near ■ According 
to his calculation, the double-core vortex is stable below 
about 2.5 MPa over the A-phase-core vortex. Combining 
his result with the phase diagram of Fig. 1 . wc naturally 
anticipate another phase transition below 0.01^02 from 
the A-phase-core lattice to the double-core lattice at low 
pressures. 

To find the transition, I have performed a variational 
calculation down to Cl = 0.0001Oc2 at p = 0.0 MPa us- 
ing N < 3000 Landau levels. The hexagonal lattice is 
assumed to simplify the calculation, and I have taken 
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FIG. 8: Order parameters along x axis at p = 0.0 MPa and 
n = 0.00015nc2. (a) O < a;/^ < 10. (b) < a;/^ < 100. The 
amplitudes arc normalized by the f2 = value of A^x, and 
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0. 



FIG. 9: Order parameters at p = 0.0 MPa and Sl = 0.0001f2c2. 
(a) Along X axis over < a:/^ < 10; (b) along y axis over 
< < 10. The amplitudes are normalized by the f2 = 
value of A^x, and the common phase 7r/8 is subtracted from 
Ay,i. Notice ReA^x(0)^Iml2;;(0) and ReAa,^(0) ^ImAL;,^(0), 
indicating the double-core vortex. 



72 X 72 integration points of equal interval per unit cell. 
It has been checked that increasing integration points 
further does not change F beyond 10^^ order. Unfor- 
tunately, however, N < 3000 Landau levels at 0.000inc2 
are still not enough to obtain enough convergence. In- 
deed, it has been observed that including more Landau 
levels as TV < igQO, 2400, 3000 change the amplitudes 
of the core superfluid components to a noticeable level 
and lower the critical angular velocity of the transition. 
Also, since the isolated double-core vortex has only a 
discrete 7r-rotational symmetry about z axis, the double- 
core lattice is expected to deform spontaneously from the 
hexagonal structure. Thus, the present calculation is a 
variational one to clarify the nature of the transition as 
well as to estimate an upper bound for the critical fi. It 
should be noted, however, that assuming the hexagonal 
lattice hardly change the free energy in this low-Q region 
with a large intervortex distance, and does not affect the 
critical fl if the transition happens to be second-order. 

Figure 8 displays the order parameters along x axis 
at p ^ 0.0 MPa and n = 0.00015f7c2. The relations 
Rei^a,(0) = Imi^j^(O) and Rei^^(O) = Imij^^(O) still 
hold; hence the system remains in Phase V of the A- 
phase-core lattice. Comparing Fig. 8(a) with Fig. 7, we 
observe that the core superfluid components have grown 



substantially, and the bulk components A^x, Ayy, and 
Azz are somewhat depleted in the core region. However, 
Fig. 8(b) shown over the longer distance < x/^ < 100 
indicates that the requirement ReAxz{0) = Im^j^^(O) 
of the A-phase-core vortex is unfavorable for a further 
growth of KeAxz- Thus, we see clearly that the system 
wants to break the symmetry of the A-phase-core vortex 
around fi~0.00015Oc2. 

Figure 9 shows the order parameters for a slightly 
lower n = 0.0001ilc2 at p = 0.0 MPa along (a) x axis 
and (b) y axis. Here the symmetry of the A-phase-core 
vortex is manifestly broken as BieAxz (0) ^ ImAyz (0) and 
BjeAzx{0)^lm.Azy{0). Thus, the double-core lattice, also 
called Phase VI, is realized. Indeed, these have the 
same qualitative features as those of the isolated double- 
core vortex calculated by Thuneberg.^^ The order param- 
eters are given explicitly by 



A<fi^=VvJ24li^N^., (41a) 
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JV 

+ E[4q'>iVq, +4+4q,^iV+4qJ , (41b) 

TV 
AT 

-i\fV X^I'^A'qiWqi -C?W~+4qiV'iV+4qi] , (41c) 
JV 

4'^ = ^ Et^i^qiV'iVqx +C^l2q,V'iV+2qJ 

AT 

+ Et ^IV-iVq, +4+2q,V'iV+2qJ , (41cl) 

JV 

iW = -^^/^7 ^[c^qi'^^qi ~'=^+2qi^^+2qi] 

AT 

-C;Ar+2qiV'iV+2qJ , (41e) 

N 

= Et^i^+lqiV'iV+lq, +cj^Uq,V'JV+5qJ 

AT 

+ VF Et 4+lq,V'iV+lqi +4+5q.V'AJ+5qi] , (41f) 
N 

= isfv E[c^+iqi'^^+iqi-c^+5qiV'Jv+5qi] 

Af 

-iV^ E[4+lqi^iV+lqi-4+5qi^JV+5qi] , (41g) 
AJ 

^i"'^ = ^ E[4"+klV'iV+5q, +4"i,VAJ+lqJ , (41h) 

N 

^i'^ = ^ El^+lqi^^+l'll +4+5qiV'Ar+5q,] , (41i) 
AT 

with N = 6n, /? = 7r/3, and p=l. All the coefficients are 
real, and terms with which are absent in the A- 

phase-core vortex, bring new asymmetry between 
and Ay™-* as well as new terms in Ai'^^K 

It follows from Figs. 8 and 9 that a phase transition oc- 
cms between = 0.00015fic2 and Q = 0.0001ilc2 at p = 0.0 
MPa. As already noted, however, N < 3000 Landau levels 
are still not enough to identify the critical value q'^^^^^ 
or the order of the transition. Including more Landau 
levels has been observed to decrease ric^^^^-*, so that 
the present value fi^^"^^' « O.OOOir2c2 should be con- 
sidered as an upper bound for oi^""^'^ at p = 0.0 MPa. 
It is also expected from Fig. 1 that is a de- 

creasing function of p. As for the order of the transition, 
the isolated- vortex calculation by Thunebcrg^^ clarified 
that it is first-order near Qci- However, it can be second- 
order from a purely group-theoretical viewpoint. There- 
fore, it is possible that the transition at p = 0.0 MPa is 
second-order, changing its nature into first-order through 
a tricritical point '^'^ as we increase p. 

Experimentally, the present result implies that, as 
we increase n above 3 rad/s, the double-core region 
will shrink in the p-T plane to be replaced by the A- 



phase-core region, vanishing eventually in a rotation of 
10^ ~ 10^ rad/s. In addition, the phase boundary be- 
tween the double-core and A-phase-core vortices may 
change its character from first-order to second-order in 
increasing O. A systematic study for f^<10^ rad/s may 
be able to detect this shrinkage of the double-core region 
as well as to provide an estimate for the critical Q where 
the double-core vortex disappears. 



VII. SUMMARY 

Based on a new approach starting from flc2 down to- 
wards flci, the present calculation has revealed a rich 
phase diagram of rapidly rotating superfluid ^He. Six 
phases have been found in the p-Q plane, and we now 
have a complete story of how the polar hexagonal lattice 
near Oc2 develops into the A-phasc-corc and double-core 
vortices experimentally observed in the B phase. Inter- 
estingly, the polar or the ABM state is favorable over the 
isotropic Balian-Werthamer state for f7>0.ir2c2- 

The present study have also clarified prototypes of vor- 
tices expected in multicomponent superfluids and super- 
conductors. With multiple order parameters, it is pos- 
sible, and may be energetically favorable, to fill cores of 
an existing component with others not used yet. These 
unconventional vortices have been classified here into 
two categories, i.e., "shift-core" and "fill-core" states. 
The former vortex lattices are composed of interpen- 
etrating sublattices of different components which can 
be described by using multiple magnetic Bloch vectors 
(qi, q2, • • • ). They have an enlarged unit cell with mul- 
tiple circulation quanta. This category includes Phase 
III of Fig. 1 where the mixed-twist lattice is realized 
at low f2, and also, the vortex sheet expected in two- 
component systems. In both cases, spatial variation 
of the 1- vector is the main origin of vorticity. The latter 
vortex lattices, which may be realized for 0<fic2/3, can 
be described using odd N Landau levels with the same 
magnetic Bloch vector. This includes the A-phase-core 
and double-core lattices in the B phase. 

Superfluid ^He is a unique system with 9 order pa- 
rameters without intrinsic anisotropics. In spite of every 
difficulty to realize r2c2~ (1— T'/Tc)xlO'', the whole phase 
diagram over < < ilc2 is worth establishing exper- 
imentally to advance our knowledge on unconventional 
vortices. As a first step of this project, however, it may 
not be so difficult to observe shrinkage of the double-core 
region in the p-T plane in increasing Q. Observations of 
vanishing double-core region in the p-T plane and the 
appearance of the normal-core vortex will mark next two 
stages. To realize f2c2, it will be necessary either to ac- 
quire n ~ 10^-10^ rad/s at low temperatures or to per- 
form accurate experiments in the region 1—T/Tc^l0~^- 
10~^. In the former case, the sample size should be made 
adequately small to keep the pressure constant over it; for 
example, Ap = 0.1 MPa at = 10^ rad/s for the sample 
radius of i? = 5x 10~^ m. 
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Theoretically, it still remains as a future problem to es- 
tablish a complete phase diagram of the A-phase region, 
i.e., how the mixed- twist lattice of Phase III grows into 
either of the locked vortex 1, the continuous unlocked 
vortex, the singular vortex, or the vortex sheet, which 
have been established by low-Jl calculations.^ This fact 
also implies that there may be multiple unknown phases 
waiting to be discovered experimentally in moderate fl. 
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APPENDIX A: PROPERTIES OF ^^A^q 

(i) The basis functions satisfy the orthonormality: 

(V'JVqIV'JV'q') = ^JViV^qq' 



X]l^W(V'JVq| = 1 



(Al) 



(ii) Upon applying Eq. (18) with R = riiai +n2a2, the 
basis function i^Nq is transformed as 



TnipNqir) = e 



— iq-R— i-n-nin2 



V'JVq(r) . 



(A2) 



(iii) The function V-'Wq is obtained from ^jvo by a mag- 
netic translation as 



V'JVq(r) = Ti2 xitpNo{r) ■ 



(A3) 



It thus follows that tpNq and ipNq' are essentially the 
same function, differing only in the location of zeros, 
(iv) The basis function ■^jvq(r) vanishes at 

R rnin2=odd; A^: even, 

r=^-^cqx=2 for <^ (A4) 
^ I nin2=even; N: odd. 



Thus, ■i/'Arqi(O)=0 for qi = (bi-hb2)/2 and A' = 2n. 
(v) Generally, VArqi(r) and ipNo{^) at least satisfy 



ipNqi{~r) = e 



V'ivo(-r) =e'^"Vivo(r) 



(A5) 



(vi) For centered rectangular lattices with p= 1, ■0Arq(r)'s 
for q = and qi satisfy, in addition to Eq. (A5), the 
following equality corresponding to the mirror reflection 
with respect to the plane including z and ai-|-a2: 

ipNq{—xcosl3 + y sin j3,x sin P + y cos/3) 

= exp{-i(pNq) 'ip*Nq{x, y) , (A6) 
where <^jvq is a constant which does not depend on r. 

(vii) For the square lattice with (3 = tt/2 and p= 1, the 
basis functions tpNq^{r) and tpNoi^) satisfy 



f VjVqi(E;'r) =6^(^-1)'^ Vjvqi(r) 
V'JVqi {x, -y) = e''^/2 V^q^ (a;, y) 
Vivq, (y, a;) = e--(^-2)/2 (a;, y) 
V'JVo(E;'r) = e^^'^Vivo(r) 
ipNo{x,-y) = tp%o{x,y) 

[ Vjvo(y,a;) = e--^/2V]Vo(a=,2/) 



(A7) 



where denotes a rotation around z axis by ip = m:/2. 
It hence follows that V'Afqi(O) = except N = 4:n+1, and 
V'Jvo(0) = except N = 4n. 

(viii) For the hexagonal lattice with f3 = n/3 and p= 1, 
the basis function ^^^^{r) satisfies 

' Vjvqi(E;V)=e'(^-i)*'Vjvq,(r) 
< V'JVqi(x,-y) =eW4^^^^(a;,y) ^ (A8) 

^ Vivq, {-X, y) = e-/4+'(^-i)- rNq, i^, y) 
where ip = n'K/3. Thus, V'JVqi(O) = except A'' = 6n+1. 
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